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Abstract 

We study general gauge-dependent dynamical equations describing homogeneous 
isotropic cosmologies coupled to a scalar field (scalaron). For flat cosmologies (fc = 0), 
we analyze in detail the previously proposed gauge-independent equation describing the 
differential, x(a) = of the map of the metric a to the scalaron field ip, which 

is the main mathematical characteristic locally defining a ‘portrait’ of a cosmology 
in the so-called ‘a-version’. In a more habitual ‘p-version’, the similar equation for 
the differential of the inverse map, x(b) = X~^(q^)) can be solved asymptotically or 
for some ‘integrable’ scalaron potentials v{ip). In the flat case, x{ip) and x(®) satisfy 
the first-order differential equations depending only on the logarithmic derivative of 
the potential, l{ip) = v'{ip)/v{ip) . Once we know a general analytic solution for one of 
these x-functions, we can explicitly derive all characteristics of the cosmological model. 

In the a-version, the whole dynamical system is integrable for k ^ 0 and with any 
‘a-potentiah, v{a) = v[ip{a)], replacing v{ip). There is no a priori relation between the 
two potentials before deriving x(®) or x{ip), which implicitly depend on the potential 
itself, but relations between the two pictures can be found by asymptotic expansions or 
by inflationary perturbation theory. We also consider alternative proposals - to specify 
a particular cosmology by guessing one of its portraits and then finding (reconstructing) 
the corresponding potential from the solutions of the dynamical equations. 

The main subject of this paper is the mathematical structure of isotropic cosmolo¬ 
gies, but some explicit applications of the results to a more rigorous treatment of the 
chaotic inflation models and to their comparison with the ekpyrotic-bouncing ones are 
outlined in the frame of our ‘a-formulation’ of isotropic scalaron cosmologies. In par¬ 
ticular, we establish an inflationary perturbation expansion for x(®)- When all the 
conditions for inflation are satisfied, which are: u > 0, A; = 0, < 6, and x(®) 

obeys a certain boundary (initial) condition at a ^ —oo, the expansion is invariant 
under scaling of v and its first terms give the standard inflationary parameters, with 
higher-order corrections. When v < 0 and 6x^ < 1 our general approach can be applied 
to studies of more complex ekpyrotic solutions alternative to inflationary ones. 


* Alexandre. Filippov@jinr. ru 
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1 Introduction 

In this paper we mostly study general Friedman-type isotropic cosmologies with one effec¬ 
tively massive scalar held {scalaron, or, inflaton), see, e.g., |I] - |3]. First of all, we have in 
mind cosmological models pretending to describe a ‘pre-Big-Bang’ evolution of the Universe 
in the frame of general relativity supplemented with a scalar held the exact nature of which 
is not important for us. The most successful models of this sort were discovered in papers 
H-i in which there were proposed various models of inhationary cosmologyli] However, at 
the moment, some alternative models are not completely excluded by the observational data 
and also attract attention (see, e.g., papers na, na, im - [19] , and many references therein). 

good supplement to its general presentations in books m - 0 might be the summary of main ideas 
in m and of more recent development in m-m- For our notation see also Appendix 6.1. 
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Here our aim is to spell out the mathematical structure of the general isotropic cosmolo- 
gie^ and to study, in particular, gauge independent solutions of the cosmological dynamical 
systems in the a,^jJ,t-versions. The last one, which we call the standard version, is usually 
considered as a most useful in cosmological consideration. However, the dynamical equations 
are rather complex even in the case of vanishing curvature parameter k, when they can be 
analytically solved only for simplest potentials. For this reason, there was invented (after 
many ‘trials and errors’) an intuitive ‘slow-roll’ approximation which allowed one to con¬ 
struct the so-called chaotic inflation, describing a very short but important interval in the 
evolution of homogeneous isotropic cosmologies [7]. This approach to cosmological evolution 
is rather visual and agrees with all well established cosmological data but it leaves aside the 
initial conditions and global picture of classical solutions, not speaking on existing viable 
alternatives to inflation. In our former search for exact integrals and approximations in 
non-isotropic cosmological and static reductions of spherically symmetric gravity coupled to 
scalars (see [20]-[29]), we had found a different formulation of the relevant dynamical systems 
which we here develop and apply to general isotropic cosmologies including inflation. 

The hrst nontrivial example was the integrable model of gravity coupled to a massless 
scalaron studied in [20]. The idea is to construct the potentials, for which there exist addi¬ 
tional integrals. Then the a-dependence of the dilaton function f, which is the part of the 
higher-dimensional metric, can be derived by one integration. For the simplest model one 
can draw the (a,/9)-portrait of the static and cosmological solutionsjfl Like a phase portrait 
of a simple dynamical system, it consists of curves Filing a domain in the (a,/5)-plane and 
having a few singularities which look like nodal or saddle points. We derive the curves of this 
portrait analytically but in general it may be considered as an object of differential topology. 

This seemingly artihcial ’’looking for the lost keys under the street lamp” turned out 
to be useful both in constructing approximate solutions [21] and in hnding new integrals of 
motion [27]-[29]. The main problem in anisotropic cosmology and in general static models is 
that we have three dynamical equations, two for the metric functions a{t) and f3{t) and one 
for the scalar matter ff), while we are granted only one integral - the energy constraint. 
A symmetry of the potential sometimes provides us with one more integral but the third 
canonically commuting integral is a rarity, at least, in physically interesting systems. 

The structure of the spherically symmetric reduction, which is the two-dimensional held 
theory, allows one to hnd some integrable classes of models if we make strong simplifying 
assumptions about their potentials. For some multi-exponential potentials and the simplest 
(‘minimal’) coupling of scalars to gravity, there exist intemable systems related to Liouville 
or Toda-Liouville two-dimensional theories (see H-1251)1 □ For the one-dimensional cosmo¬ 
logical reductions with one scalaron there might exist more integrable models. The simple 
exponential potential and a special potential with two exponential terms are integrable. But 
more interesting polynomial potentials are considered as not integrable and are usually stud¬ 
ied in approximations that can be treated asymptotically or numerically. What must be the 
scalaron potential in the models of very early Universe is not completely understood. In the 
‘standard’ model of chaotic inhation, the potential v{'ip) proportional to 'ip'^ was most popular 

^We do not suppose that the space curvature parameter k is zero, do not choose a special frame (gauge), 
but for the moment assume the minimal scalaron coupling to gravity and arbitrary potential v{tp). Non- 
minimal cases can be considered by applying the well-known Weyl transformation of the metric. 

^In this model the effective potential includes the electromagnetic term. When the scalaron coupling is 
minimal, the result can be written in terms of elementary functions. A topologically correct sketch of the 
portrait derived in m (see also m) can be found in unpublished report m- 

“^Such integrable models were, in particular, obtained in supergravity theories (see, e.g., [3l]-[35]). Pure 
exponential potentials are interesting theoretically but were not so popular in cosmological applications. 
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but it may be replaced by more complex potentials. Strictly speaking, the models of the 
very early Universe, including the most successful inflationary models, do not yet form a 
completed theory, and they deserve deeper studies together with their alternatives. 

As a hrst step we propose to develop a fresh view of the mathematical structure of a 
fairly general homogeneous isotropic cosmologies. The next step will be to analyze them 
together with the (non-isotropic) static and wave sectors. The third step must be to study 
their interrelations in the frame of the two-dimensional theory of spherical gravity coupled 
to scalarons. This program had been attempted to partly implement in [22]-[27], in the con¬ 
text of the multi-Liouville and Toda-Liouville integrable models, where the three sectors of 
the solutions were derived and classihed. The difficult problem that was not really touched 
in that papers is the derivation and study of the solutions in presence of small perturba¬ 
tions. Here we make only the first step - formulating different approaches to constructing 
approximate and exact solutions of cosmological dynamics with arbitrary potentials 

The main instrument is the gauge-independent x-eguation that was hrst introduced in |36] 
for the case k = 0. One of the main results of this paper is the derivation of some exact and 
approximate solutions of the equation for x{'ijj) = afif) and a fairly complete investigation 
of their asymptotic properties for large and small if within a physically important class of 
the potentials. To gain access to a most general picture we hrst introduce and discuss the 
dynamical equations in the arbitrary linear gauge (c-gauge) dehned by the simplest choice of 
the Lagrange multiplier e'*' in terms of the ‘physical’ isotropic metric e^“, namely, 7 = —ca. 
Diherent gauge choices are hrst illustrated by integrable models and special solutions on 
which we also illustrate the concept of portraits in cosmology, generalizing the phase portraits 
for the gravitational and matter subsystems and those proposed in Refs. ra. m- 

The most informative portraits in cosmology are the differentiable maps ifi^a) and a{if) 
that can be derived by integrating x(a) = dip/da and xi'f) = da/dip, respectively. The 
closed hrst-order diherential equations for these functions can be written only when fc = 0 . 
But even then we can analytically solve the equation for xi'f) either with a few integrable 
potentials v{ip) or in the form of power-series or asymptotic expansions. In contrast, we 
show that the eguation for x^(a) can be explicitly solved if we formally replace v{ip) by 
v{a) = v[ip{a)] which would be really possible if we knew x(a) = ip'{a). Of course, we don’t 
know it before solving the y-equation. Therefore, we should act as follows: take an arbitrary 
‘potential’ v{a), hnd the exact expression for x(q:), derive its integral ip{a), take the inverse 
map a{ip), and thus eventually get v{ip) = v[a{ip)]. This is not a transformation of a given 
function v{a) into a unique function v[ip{a)] or vice versa. More precisely, in this chain a{ip) 
depends on two arbitrary integration parameters - one is an initial or boundary condition 
Co for x(a;) (see Section 4.3), the second, ipo emerges in integrating x(q;) = ip\a). In general 
it looks like a sort of ‘anchoring’ potentials in one version to some particular solutions in 
another one, i.e. the potential v{ip) will depend on Co (presumably, we can neglect ipo). @ 
Luckily, we live in a Unique Universe (the best or the worst) and Co must be fixed. In 
our general a-version formulation of the inflationary cosmology in Section 4.3 the condition 
Co = 0 gives the unique solution y using which we freely jump from the exact a-expressions 
to approximate (perturbative) ones in the ■^-version, with arbitrary potentials v{a) or v{ip). 
In this way we derive the general inflationary solutions. We also discuss a more radical 
proposal: the potential is not a fundamental input in cosmology and we may venture instead 
to using one of the portraits. From the theoretical point of view, one of the best candidates 
to this role are the gravitational phase portrait a{a) and scalaron a-portrait ip{a). 

®See the simple examples and exact solutions in Appendices and in the main text. 


4 



2 Dynamical equations 

In this paper, we consider the reduction of the two-dimensional theory to equations describing 
isotropic cosmology and ignore other one-dimensional reductions studied in our previous 
work [2D] - [25]. The procedure and notation is briefly described in Appendix 6.1. Note that 
the effective potential and the equations can also be directly derived from the spherically 
symmetric Einstein equations using the Lagrangian 

i'"’ = •/=§ - vW - . (1) 

2.1 Gauges and gauge independence 

Supposing that the matter held depends only on time t and applying the reduction procedure 
we hnd the effective cosmological Lagrangian (see Appendix 6.1) 

£( 2 ) _ e3a-7(^2 _ _ g3a+7^(^) _ Q^ga+7 ^ ( 2 ) 

where now a ,/?,7 depend on t and (3 = a. We see that here e'^ is the Lagrange multiplier 
related the parametrization invariance of the Lagrangian. We hx the class of possible gauges 
by the linear condition 7 -(- ca = 0 , where c is an arbitrary real constant. 

This class of gauges includes: the standard (ftime'’) gauge (S) c = 0, 7 = 0; the Hamil¬ 
tonian gauge (H) c = —3, 7 = 3a; the light-cone gauge (LC) c = —1, 7 = a. Here we hrst 
write and try to solve the main equations in the general linear gauge, 7 = —ca: 

C, = - bd") - . (3) 

First we vary Cc in the gauge parameter 7 and write the Hamiltonian (energy) constraint: 

Hc^ 7 ^ — 6 -I- ^(-0) + Qk = 0; r] = ip , ^ = a , (4) 

where we introduced the momentum-like variables 7 , ^ that are more convenient than the 
canonical momenta = dr^Cc and pa = d^Cc ■ Note that the canonical Hamiltonian, 

^can _ ^-(3+c)a|py4 ^ ^ g Pa = 12 ^6^^+ (5) 

coincides with the constraint He only in the gauge c = —3, and we see that H'f'^ can never 
be split into dynamically independent scalaron and gravity parts (unlike Hf). 

The first-order eguations (equivalent to the canonical ones) are derived by varying La¬ 
grangian (|3]) in fj, a, then replacing ip, a by 7 , and, hnally, using constraint (0]): 

Ip = 7 ], d = G 2 r] -\- 2{3-\-c) T]P, v'(ip) = 0 , ( 6 ) 

6 ^ + (3 + c) 7 ^ + c e-2=“ v{ip) + 6k {1 + c) = 0 . (7) 

Expressing the n-term in the last equation using constraint (jlj) we can rewrite it as 

2 ^ + 7 ^+ 2 ce^+ = 0 , ( 8 ) 

which is independent of v{ip) in any gaug^and gives several interesting exact relations. For 
instance, taking in it c = /c = 0, we obtain below a very simple and useful equation fl2l|) . If 

®In the S-gauge, c = 0, it coincides with the previous equation ©• 
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c,k>0, we find the general exact ineqnality ^ = d < 0 . Similarly, from Eq.([n]) in the gange 
c = —3 we hnd ?) = -0 < 0 if n'('0) > 0. One can derive from (|8]) other interesting ineqnalities 
and exact relations for different valnes of c, k and independent of Moreover, one may 

consider equations (jl]), ([H]) as the fundamental system completely dehning the cosmological 
solutions. The evident application is the reduction of the whole system to one differential 
equation for ^(-0) in the case k = 0, see (j38il . Below we also propose a much less evident 
approach based on reinterpreting equations (EJ-dH]) as equations for ^^(a), rf{(y). 

These equations form our dynamical system of the hrst-order differential equations while 
'Hci'di '0! *^) is their Erst integral, which is constrained to zero as required by the parametriza- 

tion invariance. The vanishing of "He is a separate condition if we wish to forget the origin of 
equations (I6])-(IH1) from the Lagrangian theory ([ 2 ]). The standard formulation of cosmology 
uses the S-gauge and the second-order form of equation ([ 6 ]) supplemented by Eq.(jl]). As we 
here use all equations dl])-© in different gauges and with different parameterizations of the 
dynamical variables it is important to clearly understand their interrelation. 

If equations dSD-dZI) s^^e satished, 1-ic is their integral of motion, i.e., "He = 0 on their 
solutions. One can also hnd that, on solutions of equations dSD and dH]), the constraint 1-ic 
satishes the equation Tic = —2cf'Hc- This means that their integral of motion is = Td-, 

which coincides with the constraint if c = 0 and thus "Hq is the integral in S-gauge. Now, 
suppose that we solve constraint (jl]) and one of three equations d6])-dH])- Then it is easy to 
check that the two other equations are satished by the solutions of the chosen pair. 

2.1.1 Gauge-invariant eqnations and general remarks 

Before we pass to more concrete problems we write the above equations in the gauge- 
independent form, which happens to formally coincide with the S-gauge. Indeed, equations 
dlD-dH]) are easily transformed into the c-independent form if we dehne the gauge invariant 
momentum-like variables fj, and f and evolution parameter r. 


dr = e dt, d/dt = 

e “ d/dr , p = e ^ rj, f = e f , 

(9) 

in terms of which the main dynamical equations acquire the gauge-independent form: 


dip/dr = 7 ], 

2 dr]/dr 6 pf v'{ip) = 0 , 

( 10 ) 

da/dr = f, 

2dl/dT + f + 2 fce- 2 “ = 0 . 

( 11 ) 

H = e^^^Hc = 

ff — + v{ip) -I- = 0 . 

( 12 ) 


These equations are identical to Eqs.dlj), (El), (0) f^e gauge c = 0 , if we omit bars and 
identify d/dr with d/dt (i.e., with the dot diherentiation). In what follows we use this 
interpretation of flTOll - flT^ without comments, if this will not lead to misunderstanding. 

Here we consider the simplest equations for spherically symmetric gravity coupled to a 
scalar held. As is well known, the simplihcation is dehned by the isotropy condition (see 
fl9^ in Appendix 6.1), which requires the equality (3{t) = a{t) that is in fact an additional 
integral of the complete dynamical equations. To solve the reduced equations (ED-dZI) we need 
just one additional integral that commutes with the constraint. For instance, if = 0 we 
see that, in the H-gauge, rj is constant and the equations of motion can be explicitly solved. 
In anisotropic cosmologies, there is an additional second-order equation for (dit) and we need 
one more integral to solve the equations of motion, see the examples in [20], |27]-|29]. 

The standard approach to cosmology usually deals with equations for time-dependent 
functions like (E|)-(I7j), or, with their second-order form. Also, the standard treatments uses 
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the gauge c = 0 ignoring all possible gauge transformations. Here, we mostly consider the 
general gauge non-invariant equations with ip or a treated as independent variables. As 
mentioned in Introduction, to solve the theory it is sufficient to derive either xii^) = 
or x(a) = ip'{a). We show that, for fc = 0, it is possible to hnd a hrst-order differential 
equation expressing x'(V’) in terms of a third-order polynomial of y with coefficients that are 
rational functions of ip if, for instance, v{ip) = e^°‘R{'ip), where Riip) is a rational function. 
The similar equation for x(a) can be completely solved in an unusual way that is briefly 
explained in Introduction and will be explicitly demonstrated below. To compare these two 
approaches we also discuss some exact and general asymptotic solutions of the x equation. 

The characteristic features of our approach are the following. 1. Unlike the standard 
practice, we use different gauges and versions of dynamical equations, in which dynamical 
variables are parameterized not only by the gauge dependent parameter t, but also by in¬ 
variant (‘physical’) parameters ip or a. 2. In addition to the standard phase portraits, ri{ip) 
and ^{a), we study ‘twisted’ ones, 7]{a) and ^{ip). 3. The most important for us are the 
locally equivalent, invariant maps (portraits) ip{a) and a{ip) as well as their differentials 
(‘sketches’) x(a) = dip/da and y = da/dip. 4. We do not spent much efforts on search for 
classically integrable potentials; instead we try to hnd approaches to constructing exact and 
approximate solutions describing the portraits of wide enough classes of cosmologies. 

To illustrate these features we begin with a brief review of known integrable examples. 

2.2 Simple examples from ‘upside-down’ standpoint 

With different forms of the equation and using different gauges one can easily solve several 
special cases. The most obvious integrable cases are: the constant (‘cosmological’) potential 
V = Vq = 2A and the more complex exponential potential, v = Vq exp gip. They are important 
in the context of our approach and will be considered from different viewpoints. In the hrst 
case, the last equation in ([ 6 ]) is easily solved in any gauge: noting that g/f = dp/da we 
hnd T] = r 7 oexp[—(3 -f c)a]. In the H-gauge, this integral of motion is simply rj = tjq; in the 
LC-gauge we have the integral e^°‘r] = t]q, which corresponds to the well-known integral of 
motion (fip = = Cq of non-isotropic (/9 ^ a) models of cosmologies and static states 

with the -^-independent potentials0 Now, having the explicit expression for p{a) we can 
derive from the constraint (| 1 ]) and thus hnd tc = f da/a(a) with a(a) dehned by: 

= Vo -h 6k -|- pg ■ (13) 

This is the gravitational phase portrait d(Q;) of the cosmology with the constant cosmological 
potential. To derive the ip{a) portrait is easiest in the H-gauge: when c = —3, we have 
ip = Po (t — to) and thus the gauge-independent portrait ip{a) will be found if we derive t{a) 
in the same gauge. The explicit general expression for t{a) can be obtained from fll3p but 
for simplicity we write the hnal result for A: = 0 , -Uq > 0 (for Uq < 0 , sinh^ i—)■ cosh^): 

77oe"®“ = Vo sinh^[vW2ho {t - Ao)] = Vg smh‘^[^/^{ip - ipg)] . (14) 

Like the portraits of simple integrable dilaton gravity models discussed in [20], [27], this 
portrait essentially depends on one free dimensionless integration parameter Pg/vg but has 
no interesting singularities except the usual ones ata^±oo,'^^cx),'^^' 0 o- 

^The important property of this integrable dilaton gravity model is that for Cg 0 there is no horizon, 
which reappears when Cg = 0. If the potential depends on the dilaton field y, it is in general not integrable. 
For integrable potentials n(y>) we derived two-dimensional portraits that are the systems of curves in the 
plane [h = y), which essentially depend on one parameter Co , and look like a phase portrait ([20], [27]). 
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Naturally, the model with the constant cosmological potential can be simply solved and 
analyzed in any gauge. The search for integrability of more complex cosmologies can be 
simplihed by other gauge choices. Though the S-gauge is preferred by cosmologists, the 
H and LC gauges are often more convenient. The LC-gauge is indispensable in studies of 
unihed description of static and cosmological models and it most directly relates them to 
higher dimensional theories (see, e.g., |2n]-[22], [26]-[29] and references therein). The H-gauge 
was extensively used in our searches for integrals and integrability of these models. 

2.2.1 Solutions of equations with exponential potentials 

Especially interesting and simple example of application of the H-gauge, c = —3, is the A; = 0 
case. Then the Lagrangian and second-order equations of motion are 

Cc = - 6a^ - ; 2 ^/) + = 0 , 2 d - e®" n(V’) = 0 . (15) 

If n = Uoe®^, there exists the obvious integral of motion 'll) + g a = Cq and, as is well known, 
the system is integrable. To explicitly integrate it we dehne (p = g'lp + 6a, which satishes the 
Liouville equation having the integral -|- 2hoe‘^ = Cf, where ho = Vo{g^ — 6)/2. It follows 
that the complete solution for the exponential potential can be written as: 

^ ga = Co{t - to), + = hoCr^{l+ ecosh[Ci(t-ti)]} ; e = VQl\vQ\. (16) 

Eliminating t we hnd the relation between a and V’ that implicitly dehne the portrait 'ip{a)\ 

g-(5h + 6a) ^2hoC'f2cosh2{C'i[V' + (7a + C'o(to-ti)]/2C'o}, £ =-|-1. (17) 

It is easier to hnd from flTB]) the explicit expressions for and a{t). Deriving the inverse 
of a{t) we can hnd the ‘explicit’ representation for 'ip{a), but it seems that in this simple case 
the ‘parameterized portrait’ {a{t), ipit)) is more convenient for cosmological applications]^ 
A more interesting and realistic model described by Eqs. ffTSj) is proposed in [29]. We 
remind that the potential v = vi + V 2 is integrable when Vi are arbitrary and 
9192 = 6. Indeed, dehning new helds by the pseudo-orthogonal transformation 

\/6 a = ctjji — S‘ip 2 , 4’ = ~s4i + c '02 ; c = cosh 6 , s = sinh 6 , (18) 

we hnd that the kinetic Lagrangian is —ipf + if 91 = V^tanh^, g2 = 6/91, then 

6 q ; -|- 9 i 4 = gi4i ^ 9i = cosh 6 , 92 = \/6/ sinh 6 . (19) 

The Lagrangian in flTHll then describes two explicitly solvable Liouville models: 

C^ =+ V2e~^^^^. (20) 

As the parameters vi,V 2 are arbitrary, gig 2 > 0, and 6 may be negative or positive, it is 
possible to construct cosmologically interesting models by choosing for them proper values. 
In particular, the potential v{4) may qualitatively resemble inhationary potentials discussed 
in some models if we take Vi = —V 2 > 0 and cleverly choose 9 in the interval (— cxo, -|-cx)). 

The formal solution of this model is simple: we know 4i{t) and thus easily hnd a{t) 
and 4{t)- Then, eliminating t we can derive the a{4) portrait of the system. This looks 
very simple but details are in fact cumbersome because, in general, Si and Ci dehning the 
solutions 4iit) (by two relations like the second equation in flTbll ) may be positive or negative 
and the inverse functions tippi) are involved. All these ‘subtleties’ are essential in cosmological 
applications. For these reasons we think that this ‘bi-Liouville’ model undoubtedly deserves 
a detailed investigation that lies completely outside the scope of this paper@ 


®The explicit representations for bln) Qf('b) can easily be derived in the asymptotic region t —>■ 00 . 
®More complex exponential integrable models with one and more scalar fields can be found in [22]-[25]. 



2.2.2 Note on independence from potentials 


The above simple examples demonstrate that it may be useful to replace the time variable 
by the more ‘physical’ variable a. We also can (and will) use the scalaron variable ip, but a 
plays a special role because there exists equation ([H]) independent of the potential 
This equation can be rewritten in a very compact form: 

f] Z — — 

+ f = 0- z = ^ = = ( 21 ) 

da 

where —6z is the gravitational part of the gauge invariant Hamiltonian constraint, and 
(r)^ + v) is its scalaron part. The solution of equation (!?!]) expresses in terms of ffV^\ 



ke + z{a+) + 



( 22 ) 


where z{a+) is an arbitrary constant, —oo < z{a+) < oo. This is a nontrivial representation 
for the gravitational kinetic energy, d^, in terms of the scalaron kinetic energy, It may 
be useful for analyzing the gravitational phase portrait a{a). More generally, one may find 
entertaining the idea to use //^(a) instead of the potential v'{'ip). This resembles the simplest 
canonical transformation of p into —q and will be discussed in some detail below. 

Possibly, a better idea is to take as an input p,{a) = a{a) and derive f]{a) = 'ip{a) 
from Eq. (|Mil . In cosmological considerations, d(r) is called the Hubble parameter and its 
dependence on a is more or less understood in simple cosmological models and can be to 
some extent compared with observational data. It is possible to make reasonable guesses on 
it in various models of inflation, of bouncing, etc. From the dynamical point of view it is 
the phase portrait of the gravitational subsystem which in the domain ‘before Big Bang’ is 
the main object of interest. We have no direct observational information on this object but 
analyzing extensive literature on inflation and other scenarios of the very early cosmology 
one can develop intuition of what can be most plausible guesses for a{a). In any case, if 
we knew by any means one of the functions .^(a) or p{a), we could hnd the second function 
from equation fl2T|) and then derive from constraint dH) the potential itself as a function of 
a. To ‘transform’ it to ^(-0) we must find x(q:), as discussed below. 

Although the scalaron potential is a useful and convenient thing, we know about it much 
less than, say, on the Hubble parameter (not even speaking on the existence of the scalaron 
itself). Thus attempts to avoid using vi^pj) as an input in cosmological models look quite 
natural and this can also be done in the frame of the standard version. Then one of natural 
theoretical inputs may be the scalaron (inflaton) phase portrait 'ip{'ip) = pipj) that in turn 
can be found from some simplifying ansatz defining a concrete model (see, for example, the 
derivation of the potential in the model of ‘fast-roll inflation’ in Ref. |17j). 


2.3 Equations for x(^) their main properties 

Now we suggest to forget, for some time, the brilliant original treatments of the isotropic 
cosmology and try to have a fresh look at its mathematical structure. My Erst attempt 

central idea of this paper is to look for the portrait of cosmology, ipia) or alb), in terms of comple¬ 
mentary variables tp, a. Thus it is quite natural to consider the dynamical variables as depending on one of 
these variables instead of the evolution (‘time’) parameter. In addition, we also discuss the standard phase 
portraits, ^(a) = a{a), nW — P’Wj as well as the ‘twisted’ portraits, 77(a) = '>p{a), ^{ip) = a{ip). 

^^Here and below we often omit da, dip, dx, when possible. 
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of this sort was made in [35] (see also [37|, [33]) with the aim to understand cosmology of 
the affine generalization of Einstein‘s gravity. In that paper I have actually proposed (in a 
somewhat misleading notation gauge independent equation (33) for y = el'll}/da, regarded 
as a function of V’- We return to this equation later but hrst elaborate notation and try 
to explain more precisely our plan. The principal goal of this paper is to hnd the analytic 
portrait of cosmology in terms of either i/la) or a{^l>) and thus our main objects arJ^: 

Xia) = di>/da = ip/a = r]/^, x{i>) = da/di> = a/ip = ^/r]; y(a) y(?/>) = 1. (23) 


Evidently, the last relation is identity if we put into it 'll}(a) or a{'ip). 

To derive the equations for y(Q;) and xl'ip) from ([6D-([7|) we must use the transformation 
from the independent variable r to a or -0 as well as transformations between a and 'ip that we 
can easily perform with the help of the evident identities (implicitly applied in Section 2.2): 


d d d ■ d d d _ d d , \ d 
Using the above equations we hrst derive the gauge independent relation, 


2 ^ ^ ^y-jr] 
^ dr 


2k 

(x^-6)[y+ /'(V’)] + + 3/'(V’)], l’{ip) 


v'l'ip) 
v{i’) ’ 


(24) 


(25) 


which generalizes Eq.(33) of [36] but is not a closed differential equation for y, even if we 
replace ^ dr by da or by y dip. In the hrst case, we hnd that this expression actually coincides 
with Eq.(33) of [36], when k = 0. It can be formally written as the well-dehned equation for 
y^(Q;) if we multiply it by y(Q;) and then apply the last relation of fl2Tll allowing to formally 
consider I'l'ip) as a function of a . Thus dehning x{c()dl/d'ip = dl/da = I'{a), we have 


^ -6) (y^+ /'(«)) , y = V''(a), (26) 

which can be elementary solved for any l'{a). In the second case, we have the closed Abel 
equation for xpip), which is solvable just for a few potentials. 




X = a'(V') ■ 


(27) 


Equation fl26p can be rewritten as the nonlinear second-order equation for 'ip{a) if we recall 
that l'{a) = y(a) Z'('0) and divide it by y . Similarly, fl27|) may be treated as the second- 
order equation for apip) by recalling I'pip) = y('0) l'{a). In either transformations we should 
suppose that Ipip) or T{a) is known and thus the hrst case looks preferable. We do not study 
these second-order equations and in what follows concentrate on (j26ll and fl27ll . As vpip) is 
a more familiar thing, we begin our study with investigating solutions of (1271) . 

Note that constraint (01) gives, among other things, the important inequalities for y, y: 


y^ < 6 , 6y^ >1, ifn>0,fc>0; y^>6,6y^<l, ifn<0,/c<0 . (28) 


These restrictions on v and k are often used in cosmological considerations. If they are not 
realized, the behavior of y-functions becomes much more intricate, and we usually adhere 
to them. When y^ <C 6, 6y^ S> 1, there may exist inhationary solutions, which we discuss 

this Subsection we use the standard gauge c = 0 but omit bars over -q. 
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in Section 4.3. For negative potentials, there exist viable alternative cosmological scenarios 
including a hot (ekpyrotic) compression followed by a bounce. 

In fc = 0 case the equation for x(a) (or, xi'^)) dehne the portrait (or, Q!('^)) that can 
derived by directly integrating the solution. Moreover, once we know one of these functions 
we can in principle derive the complete solution of the cosmological equations. This is 
obvious when k = 0 but can be generalized to A; 7 ^ 0 , if we somehow find the general y. 
Then, taking in Eq. (l?I|) f]'^ = we obtain the linear equation for and so get: 

g-2a+/x2(a)j ^ (29) 

expressing the Hubble function in terms of x^(q:). By changing the integration variable and 
remembering that X^(tt) da = X~^(' 0 ) dip, a'{ip) = xii^) we hnd in the -^-version: 

pp^PiP) = j . ( 30 ) 




In applications of these formulas one has to carefully dehne the limits of integrations and 
the arbitrary constants to guarantee the positivity of As a hrst application we note that 
by using fl2^ . (I5U]) . (IT^ one can derive the potentials v{a) and v{ip) in terms of y or y: 

v{a) = [6 — y^(a) ] ^^(a) — = [6 — x~‘^{ip) ] ^^("0) — = v{ip). (31) 


Also note that Eq. fl30ll can be used in iterations involving xi'P’) and or rj, see Section 3.1. 

Equations fl2^ - fin|l demonstrate (for arbitrary k) the most important property of the 
functions y, y - once we know one of them, we can derive not only the portrait {a,ip) but 
all characteristics of the cosmological solutions: r(a), riip), v{ip), v[ip{a)], gravitational and 
scalaron energies. Even more interesting for applications may be the expression for y^(a;) in 
terms of the Hubble parameter as a function of the metric, H{a) = f{a). It can be obtained 
by rewriting equation ( 1 ^ that was also used for deriving Eq. (l 2 ^ (recall that da = Hdr): 


X 


2 


fd\nH 


ke 

) 


2 

IP 


dH , 

— - \-ke 

dr 



(32) 


This simple ‘inverse’ formula allows us to reconstruct the complete solution of all dynamical 
equations and the potential in the interval of a or r in which we know H. 

It is worth noting that equations (l25ll - (l2711 have an unusual property of being dependent 
on the logarithmic derivative, l'{ip) = [lnn('^)]', not on the potential itself. This property 
is shared by the differential equations both for y(Q;) and y(V’)- This means that apparently 
small deformations of the potential v{ip) may result in signihcant deformations of the y- 
maps. For instance, v = e + ip'^, n > 1, then for -0 —)■ 0 the ‘potential’ l'{ip) vanishes if 
e ^ 0, but it is singular for e = 0. On the other hand, the potentials ip'^ give very similar 
y-maps as V = n/ip. In the standard approach, n = 2 looks preferable as it gives, in the 
S-gauge, a linear equation for ip, which however depends also on d = .^ and is not a closed 
equation for ip so far as d 7 ^ 0. Note also that the Higgs-type inflationary potentials may 
produce hnite-range singularities ~ l /('0 — ipo) that can be dealt with similarly to n/ip ones. 

A more fundamental problem is what to do with cosmologies having nontrivial non- 
integrable potentials and /c 7 ^ 0. It is evident that in this case the y-equations fl25|) are not 
simpler that the complete system (jl])-([ 8 ]). However, basing on the structure of solutions of 
simple integrable cosmologies we propose here a different approach to /c = 0 equations for 
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x(a) that allows us to find its simple analytical solution for apparently ‘general potential’ 
v[a) = v[ip{a)]. Here 'ip{a) is at first an unknown function that is to be derived after solving 
the x(a)-equation depending on v{a), i.e., 'f{a) = f dax{c(). 

Indeed, suppose that we have derived x(a) corresponding to the potential v{a) and wish 
to find the potential v{'ip) using the definition of /^a), 

, dlnvia) d'lb din nth’) . 

Once we know x('^) ^md, correspondingly 'ip{a), we can in principle derive the inverse function 
Replacing in a by a('0) we find = l'[a{'f)]/x[(^{'f’)] and determine ^(-0) up to 
a constant multiplier. At first sight, by writing v{a) = v[ip{a)] we fix this problem but then 
there remains one constant that enumerates the curves ip{a) corresponding to physically 
different solutions. It looks as if the solutions corresponding to the same "^-potential may 
correspond to a one-parameter family of a-potentials and vice versa. One can hope to better 
understand this relation by further studies of the integrable examples given in Section 2.2 
and Appendix 6.3 as well as of special solutions discussed in Appendix 6.2. 

2.3.1 On what is the solution 

Before turning to studies of approximate and exact solutions of the dynamical equations let 
us formulate what we call the solution of our problem. In the Liouville sense, this means 
that the complete solution {a{t), can be formally expressed in terms of integrals and 

derivatives of the potential, up to some functional inversions (at best, one can find explicit 
expressions for {t{a), t{'ijj)), for a('0), or for 'ijj{a)). In most favorable cases, we can find a 
complete phase portrait of the solution, say, the dependence of d on a (or, on -0) with 
calculable asymptotic behavior near singularities. This is possible in few cases. More often 
one can derive a differential equation for one function, like xi'f’)^ which expresses its first 
derivative, x'{'f), as a rational function of x ^-nd if. 

A cosmologically relevant example is the much studied Emden-Fowler type equations 
(see, e.g., iBg-mi, H). To illustrate their relation to cosmological models let us consider 
equation (1T0|) for ifi^r) derived in the c = 0 gauge, which is usually written as: 

3 aif + v\if)/2 = 0, ip = T], a = f. (34) 

If we suppose that d is a constant, d = , and v = aif"^ — bif^, where p is a rational number, 

we obtain the Emden-Fowler equation. Actually, in cosmology described by equations (ITOD - 
(HI, the only possibility to have a nontrivial solution of equation fl3T|) with the constant 
Hubble parameter d = is taking k <0 and n = Uq -|- aV’^, so that it becomes linear @ The 
general Emden-Fowler equation is usually rewritten in the form of a first-order nonlinear 
equation, e.g., ri'{if) + 3fo + v'{if)/2r] = 0, which is the phase portrait of equation fl3T|l . Even 
with simple rational potentials, it cannot be solved analytically. Nevertheless, we better 
avoid calling it non-integrable. Indeed, the behavior of their solutions for rational p was 
analyzed in great detail, including their exact asymptotic. More recently, a few significant 
results were also obtained on their classical integrability (see 021, 031). 

This example demonstrates why a general equation like z'{x) = R{x, z) with simple ratio¬ 
nal function R{x, z) may be regarded as essentially integrable: 1. The topological portrait of 

^^See equations (I105I) - (I107I) in Appendix 6.3, where we also discuss further examples of the potentials for 
which the system cnii-dii]) with ^ = Cq or t) = Co can be solved. 
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the solutions can be found as far as we can derive the zeroes and poles of R outside of which 
the solutions are locally analytic; 2. According to Hardy’s theorem |11]-|15], the possible 
asymptotic behavior of real continuous solutions is either 2; ~ or 2; ~ aa;^(ln2;)^/"', 

where n is an integer, P{x) - a polynomial. 3. It follows that in simple cases we can derive 
approximate analytic behavior of the solutions in the (x,z)-plane, including singular points. 

Unfortunately, in cosmology we need more detailed information on solutions. For in¬ 
stance, the inflationary behavior is apparently hidden in some subtle properties of the po¬ 
tential, probably, in complex 'll) or a plane. Moreover, solving the model discussed here is 
only the hrst step in looking for observable effects. Thus a sort of classical analytic integra- 
bility of homogeneous cosmological models with scalaron is highly desirable. Some relevant 
results on integrability in static and cosmological models can be found in [20]-[29]. 

3 Dynamics in ^-version 

Here we study in detail the ■^-version equations, which are important because the potential 
u('^) can often be determined by some held-theoretical model, in which it has a certain 
‘physical’ interpretation: the mass squared term of the scalaron, a ’Higgs-type’ potential, one 
of many possible potentials derived in reductions of supergravity (see [31]-[35] and references 
therein). In connection with inflation, a wider spectrum of potentials was discussed in papers 
mentioned in this paper. The present author recently discussed a more exotic origin of the 
scalaron and its potential from affine generalizations of Einstein’s gravity theory, |36j-|38j. 

3.1 Main cosmological equations 

We mostly discuss properties of x(V’) satisfying equation fj?711 for k = 0. Explicit general 
solutions of this equation can be derived only for very special potentials vlip). The simple 
examples with the exponential or bi-exponential potentials were treated above in the frame 
of the f-version using the explicitly integrable model ffT^ - fl20l) (see also Appendix 6.2). We 
can also solve them in the -^-version but the solution is a bit more cumbersome and we’ll 
only give the solution xlip) for simple exponential potential. A more complex explicitly 
solvable model is discussed in Appendix 6.3. For general potentials u('0), we derive asymp¬ 
totic approximations at V’ 00 and V’ 0. For small 'ip^ we also construct power-series 
expansions. 

The equations in the ■^-version immediately follow from (I6i)-(j7|): 

^ = 2^ + 2(3 + c)f + fA!,.(^)=0. (35) 

dip T] dip T] 

dp g-2(l-|-c)a 

Q — + {3 +c)t] + c - v{ip) + {1 + c)6k -= 0 , (36) 

dip Tj 7] 

the constraint is unchanged. The gauge-independent equations for ^ ’fj = rje^^ are: 

2 fi'pip) -I- 6 ^-t- v'pip) fj~^ = 0 , 2 ^'(ip) + f] + 2k = 0 . (37) 

The gauge-invariant constraint is given in fll2p . 

When k = 0, the constraint and Eqs. flTT]) give the closed equation for P,{ip): 

-6 ^\iP)+v{iP)=0. (38) 
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If we could find ^ {'ip) we would have the complete solution of the k = 0 cosmologies. Indeed, 
a'{ip) = xpip) = i/ff = tp = fjpip) = —2^'pip), and it follows that: 


apip) 


iptp) dj! ^ ^ ^ _ r dp) 

a/6 — vpip) J a/6.^^ — n('0) 


(39) 


Equation fl38|) is essentially equivalent to equation fl27|l for xi^ip) and, at hrst sight, it 
gives us no new information on discussed problems. However, it is worth of independent 
study as it contains the potential instead of logarithmic derivative. It can be easily solved in 
terms of the expansion in powers of ^p. To simplify notation we take x = ^J?>/2'^p, w = \/6e, 
vpip) = v{x) and write the equation and expansions of v and w. 


{w'{x)Y — w‘^{x) + h(a;) = 0 , v{x) = ^^VnX^ , w = y/ Wnx"' ( 40 ) 

0 0 

The recurrence relations for Wn cannot be solved in general and are rather cumbersome even 
for simple potentials, like h = Uq + Vix + V 2 x‘^. The hrst coefficients for this potentials are: 


wi = wl - Vo , W2 = 4:W^ ^{2woWi - ui), 6^3 = ^{wl - V 2 ) . 

We see that this expansion is inconvenient and this also signals that Eq. fl38|l is not the best 
starting point for solving even the simplest cosmologies. The perspectives with the large x 
behavior are even worse. In next Section we hnd that the W^quation is better and derive a 
good approximation for a: —)■ 00 . With this in mind, let us relate w{x) to z{x) = \/6x{P’)- 

z{x) = w{x) [w'^ — v{'lp)]~^ , w‘^{x) = v{x) z‘^{x)/{z'^ — 1). (41) 


The hrst equation in flTT|) is a convenient starting point for iterations. We just mention 
two most evident ideas for this. Taking as a zeroth approximation a physically reasonable 
wo{x) we hnd the zeroth approximation for z: zo{x) = [1 — v{x)/wl{x)]~'i. Then, applying 
(l30|) we can derive the hrst approximation ^ 1 ( 0 :) = exp[—2/ 2 :(j'^(a:)] and hence zi, etc... 
This procedure is best suited to the case n < 0 because then Zn is well dehned and \zn\ < 1. 

We only mention an alternative procedure that looks more complex but is applicable 
to positive potentials. It can be realized by writing the zeroth approximation for z{x) in 
terms of a reasonable phase portrait for the scalaron, rjoi^p)) = fjo^x). Indeed, in this case the 
constraint dehnes the zeroth approximation for zo{x) by ^o(^) = 1 + h(a;) r)(j'^(a;), and the 
iterations are produced by the inductive relations 


Wn{x) = exp[-2 / z^_^{x )\, fi^{x) = 


wl{x) 


■'n—l 


.X 


4(^) = 1 + 


v(x) 

tj^(x) 


n = 1,2, 3,... 


The advantage of these iterations, which are applicable to inhationary scenarios, is that rjlip) 
describes a phase portrait of the inhaton much discussed in studies of inhation. 


3.2 Exact and asymptotic solutions of x(t/)) - equation 

Now we return to equation fl271) that looks like a generalized Emden-Fowler equation for a 
wide class of potentials and is simpler than fl38ll . For some potentials it can be exactly solved 
and for a rather general potentials in can be solved asymptotically, for -0 —)■ 00 and ip 0. 
We rewrite it dehning the ‘potential’ u(x) and using the above notation for y; and ip\ 

dz/dx = {z'^ — l)[zu{x) + 1], u{x) = d\n y/v/dx , z = '/6x, x = \/3/2ip . (42) 
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3.2.1 Solution with v{ip) = 

Supposing that u[x) = g = g~^ ^ 0, ±1 we can derive the solution in the form 

2{g~^ — g){x — Xq) = ln( \z + \z — \z + g\~‘^) ■ 

It is easy to find exponentially good approximations for z{x) when z is large or approaching 
±1 or —g. These can be compared with the corresponding behavior of xii^) or x(a) which 
can be derived for exact solutions given in equations (lT6|) - (fT7|) . That construction used the 
additional integral g + g^ = Co obtained in the Hamiltonian gauge c = —3. It gives the 
invariant relation between g = and Then the integral and constraint 

1 + = C'o ?? e"^" , = l + hoC'o"^(l + ^x)^ exp(^V' + 6a), (43) 

together with relation ffTT)) define the complete gauge-independent solution of the model with 
the exponential potential. We can also derive x{c(), hnd the f-dependence of y and x, and 
thus study on this example the relations between all the three versions. Such a comparison 
of versions is also possible with the bi-exponential potential models (fT8|) - fl20|) and fl96|) - fll02p . 
For the general potentials, we will find exact solutions in the a-version, while in the ■^-version 
we have only asymptotic solutions to which we now turn our attention. 


3.2.2 Important transformation of xii^) and properties of /v{'ip) 

A convenient approach to asymptotic of z{x) is to introduce the transformation. 


Z[X] = 


1 — £€ 

— = -tanh^y{x ), 


where £ = sign(l — z'^) = sign(—u) 


(44) 


1 -f £e~‘^y 

Then Eq. fl42]) becomes simpler and more suited for a qualitative analysis {v{x) = v{'ip)): 

dy/dx = 1 — u{x) tanh^(?/), where u{x) = v\x)/2v{x) = I'[x). (45) 


This equation is not difficult to analyze when the logarithmic derivative of the potential 
has nice properties. Cosmo legists often suppose that v{x) is a polynomial that does not 
change sign but vanishes at a; = 0 like v{x) ~ The behavior of the corresponding 
u{x) at infinity is rather simple. Indeed, consider a polynomial potential, v{x) = Pn{x). 
Then at infinity u{x) = {N/2x)[l + ©(x"™)], with integer m satisfying 1 < m < iV. For 
u{x) = e^^^PN{x) .P^^(x) the asymptotic behavior is u{x) = -|- (iV — iV) ( 2 x)“^ -|- 0 (x“^). 

To classify the behavior of u{x) at x —)• 0 we suppose that h(x) can be represented as a 
convergent power series Poo{x) = ^u^x^: 1) If Uq 7 ^ 0, it is easy to show that 


m(x) = ^ M^x™ = —(n 4-1) Vn+i x^ [ 1 -4 


2vi 


ILm^ 


{n -4 2 )vn +2 
{n + l)u„+i 


(46) 


for n > 0, and if Ufc = 0 for 1 < /c < n; 2) If uq = 0, then for n > 0, and if ua, = 0 for 
1 < fc < n, we hnd the universal behavior for x —)■ 0: 


uix) = — 

X 


0 


UmX 


= 7 ^(^+ 1 ) [l + y^UmX" 

2x “ 


Ml = 


"^ 72+2 


{n + 1) Vn+l 


(47) 


The second formula is most general. It is also applicable to singular potentials m(x) having 
terms ~ x"”. The simplest behavior for m(x) give pure exponential and pure power poten¬ 
tials, and x”. All these patterns of asymptotic behavior of the potential m(x) are met 
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in cosmological models and we briefly discuss the most important properties of the solutions 
for large and small values of x having in mind inflation and other scenarios. 

Some potentials u{x) used in inflationary models (see recent reviews n, USD look, at 
first sight, different but actually fit into the above classification. For example, the potential 
derived in 051 for the original Starobinsky model is no(l — e which asymptotically 

gives — ~2 uo (yf and — —2vq/'iP for small "00 The same is true for the 

more general potentials discussed in [15], {gip). Similar potentials are discussed on 

connection with ‘ekpyrotic - bouncing models (see. e.g., Hi. BSD)- A typical potential 
looks as the following generalization of our bi-exponentials: 




(48) 


where integer m and real vq are positive (inflationary) or negative (ekpyrotic), ui > 0 if 
m < 0, and gi may have different signs, unlike our model fl20ll . 


3.2.3 Large ip behavior of x(^) 

Let us first consider y —)■ ±cx), —)■ 1. Using that at infinity tanh*^(j/) —>■ ±1 and supposing 
that v{x) = 'Y^VnX'^ at large a; is a polynomial of degree N (or, more generally, e^^v{x) with 
5^ < 6) we can easily derive the main terms of asymptotic behavior-. 


y 


= X =F In ^/\v{x)\ -I- Co -I- 0(e ^^), for x —)■ ±cx), 


(49) 


where cq is an integration constant and dependence on e is hidden in the exponential correc¬ 
tion to be derived in a moment. This formula shows that, for polynomial potentials, y ^ x 
and we can obtain the main terms of the asymptotic expansions of y, z for x —)■ -|-cx): 

da 


Vg — = Vg xW = z{x) = — tanh*^ [x — In ^/\v{x)\ -|- Cq ] -|- O(e . (50) 

(Jj (JJ 

Using this and higher asymptotic approximations we can And the corresponding asymptotic 
expansions for the portrait a{ip) and, by inversion, ip{a). In the simplest approximation 

Vg ^ = —[1 — v{x) + ...] = —[1 — 2Cov{ip) exp{—\/Gip) + ...], (51) 

dip 

where Cq is the redefined arbitrary constant. By integration we can derive the asymptotic 
portrait a{ip). Taking the popular ‘cosmological’ potential v{ip) = voip^ we And 

VGa(ip) = ipQ - Ip + Cq[ 1 - (fiip'^ + VGiP 1) exp{-'/Gip) + ...] , 

where ipo is the integration constant. The inverted expression defines the asymptotic of ip{a). 

Note that this approach to the asymptotic behavior at infinity in fact uses iterations with 
the zeroth approximation given by (jUj). To And the next one we write the exact equation. 


?/(x) = I/O + / m(x) [1 — tanh^|/(x)] = |/o-|-C q — 2 / dxu{x)e {1 + ee ^^) (52) 


the additional Section 4.3 we demonstrate that inflation cannot be considered in the domain ip ^ 1 
where exist powerful asymptotic expansions in the //^-version. The reason is that the inflationary value of 
| 2 :(x)| must be large while asymptotically it is close to unity. The adequate approach is the a-version. 
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where u{x) = l'{x) and yo{x) = x — l{x) is the zeroth approximation. Replacing y by ?/o in 
the r.h.s. we obtain the hrst approximation yi, in which we may expand the integrand in 
powers of and hnd the hrst exponentially small correction to j/o also independent of e: 

^ poo ^ 

yi{x) = X — l{x) + Cq — 2 / I'(x) +0(6“^^“). (53) 

J X 

Remembering the above dehnition, l{x) = In ^Jv{x), we hnd that in this approximation 
yi = X - + Co- J ^ +^0 “ y 

where the integral can be expressed in terms of elementary and special functions for a wide 
class of the potentials. The new approximation ^/6 ip = tanh^ yi is significantly better than 
fl5CT]) and can be further improved if necessary. 


3.2.4 Small ip behavior of xW 

The behavior of z{x) and y{x) near a: = 0 is more complicated because u{x) for x —?• 0 may 
be zero, constant or inhnite. In addition, when e = —1, the solutions y{x) must not vanish 
for X —)■ 0 if u{x) is singular or constant at x = 0. For this reason we hrst consider the case 
e = +1. Then for the regular positive potentials that vanish for x —?• 0 it is easy to derive a 
universal approximation that can be used for subsequent iteration^ 

px 

2/o(x) = [h(x)]"^/M [h(x)]^/^ + ..., if I/O <1; yo = l-u{x)yo. (55) 
Jo 

This solution is independent on any arbitrary parameter. This means that it may be either 
an enveloping solution or a separatrix (if it has no common points with any other solution). 

Instead of using iterations one can directly construct a power-series expansion. However, 
for qualitative analysis of solutions with realistic potentials the general approach of (l5^ 
may prove preferable. An important example is the parameter-independent solution for 
potentials (ITT)) . When 7^ 0, the solution of flT5]) with e = +1 must vanish when x ^ 0. 
Expanding it in the series y{x) = X]n=i we can hnd all a„ with n > 1 in terms of ap. 


y{x) 


X 

1 + Ms 


XUq 

2 -|- Vjq 


X^ { Mn 

X-- Ui - — - 

o ~h tig y 2 , tig 


tis 


3(1 + Usfi 


+ 0(x^ 


(56) 


Here we take into account the three coefficients a* and the third term in 03 is the contribution 
of the y^. This formula and its extensions to higher order terms are applicable to arbitrary 
parameters. In particular, if all = 0 and Ug = a is any real number, it gives y{x) for the 
pure power potential v = x^®". The two hrst terms in fl56D agree with corresponding Eq. (l55|) 
and both methods can be generalized to derive the one-parameter family of solutions near 
X = 0. They are based on splitting y{x) into ‘big’ and ‘small’ parts, y = yo{x) -\- yp{x), so 
that yp could either be expanded in a regular power series, like fl56l) . or satisfy a soluble (e.g, 
linear) diherential equation, like that one dehning yo in fl55l) . 

In deriving solution fl5^ we suppose that yo = 0 and then hnd the unique small solution 
yi for arbitrary potentials u{x). If the w-potential is regular, Ug = 0, we can construct the 
general solution y{x) = yo{x) -\- yfix) = Oq -l- Oix -|- a2X^ -|- ... by taking yo = Oq and applying 

^®This is a special solution of the linearized equation (l45|) . the general one to be discussed in a moment. 
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the addition theorem to tanh(ao + |/i). Then we either use for y{x) the linear approximation 
in Hi or directly expand y{x) in the power series in x. We hrst write the general expansion 
for tanh^(j/o + yi) in powers of yi, with notation to = tanh(j/o) and ti = tanh(?/i): 

tanh^ y{x) = tanh^(j/o + yi) = + {1 - tl^) ti (1 + h fg . (57) 

Applying this to constant yo = Qq and approximating ti = yi + 0{yf) we can easily hnd the 
linear equation for yi, the solutions of which generalizes (l55|l . Alternatively, we may use (l57ll 
to derive the expansion in powers of x. With arbitrary oq and oi = (1 — uoto), we hnd 

y{x) = ao + aiX-]^[uito + uoCQ'^ai]x'^-^{u 2 to + CQ'^[uiai + Uo{a2-toal)]}x^-0{x ^), (58) 

where to = tanh(ao), cq = cosh(ao), and 02 is the ^^-coefficient in this expansion. The 
expression for z{x) can be written with the help of the same formula (l57)l . Depending on 
the concrete values of the parameters it may sometimes be used for extrapolation of z{x) to 
asymptotic regions. It also is useful in comparing the small x behavior of z{x) for different 
potentials. Note that this solution can easily be rewritten for e = —1. 

The elementary formula fl57j) is extremely useful for deriving various approximations. For 
instance, we can take as yo asymptotic approximation and hnd yi by solving the linear 
equation obtained from fH5ll by linearizing tanh^(?/o + yi) in yi- Then we derive z{x) = 
— tanh^(?/o + yi) using Eq. fIST)) . This equation fl57)) may hnd most interesting applications 
in studies of the small x behavior of z{x), especially for singular u-potentials and for the 
e = —1 case, when the behavior of the solutions at the origin a: = 0 is more complex. 

In fact, to compensate the singularity Us/x, one must take yo = alna:-|-ao, where a = Ug. 
Then, supposing that yi = a^x"', oq = Ina and thus yo = ln(ax“), we can apply (1571) to 

solve Eq. fHSj) for both signs of £. Indeed, dehning oq = In a, we immediately see that 

tanffi(l/o -|- yi) = —(1 -f £a^x^“)“^ [1 — {yi + 0(|/^))]. (59) 

With this relation we can either expand yi{x) in the power series and solve several recurrence 
relations or, instead, solve the linearized equation for it. We only illustrate the hrst approach 
by writing a few terms for the simplest singular potential u{x) = Ug/x: 

y{x) = In(ax^) + x — ea?{l + 4x/3) x^ (60) 

In the beginning of this Subsection we met only one analytically solved equation, on 
which one can roughly test the proposed approximation methods. The simpler explicitly 
solved example in the '0-frame is presented in Appendix 6.3. On these examples one can 
check the precision of the approximate solution derived here in more detail. In this Section, 
we briehy described results of our studies in the general -^-version. This in principle would 
allow us to hnd, in main models of cosmology, the approximate {a, 'ip) portraits which are 
exact in the asymptotic regions. A more complete presentation would be possible if we 
concentrated on specihc inhationary or ekpyrotic-bouncing models not contradicting to the 
present observational data. This task requires quite diherent means and exceeds the author’s 
space-time capabilities, but in next Section we derive exact solutions that allow us to hnd 
a very simple picture of inhation. In addition, this will help us to better understand and 
derive inhationary solutions y{x). Also note that the general ‘a-solutions’ can be used in 
analyzing some ekpyrotic-bouncing cosmological scenarios. 
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4 Dynamics in a-version 

Here we present an unusual general approach to generating cosmological solutions. It is 
based on simple examples of relations between the solutions in the ■^-version and the a- 
version (see flT^ . fl22l) . (I105l) - (lll0p ). The close mathematical connection between the two 
versions is given by the functions x('0) and x(q!) introduced and discussed above. Now 
we hrst derive the exact analytic solution of equation fl26|l for with arbitrary l{a). 

Then we write the exact analytic solutions of all dynamical equations, formulate a general 
approach to inflationary phenomena and compare it to the standard one. 

4.1 Exact solution of x^(ct)-equation for k = 0 

We start by solving the k = 0 equation because to solve the general one, 

^ = (y - 6) (y+?(«))+. (et 

one has to know When fc = 0, we the obvious general solution is simply 

X^(a) = 6 — e®" h(Q!) [Co + J e^°‘v{a)~\ k = 0, (62) 

where Cq is an arbitrary constant. This solution can be verihed by inserting 
equation (|26|) and taking account of the above dehnition I (a) = Inh(a). The main problem 
with this solution is to hnd the conditions for positivity of the r.h.s. of Eq. (l6^ . We discuss 
it on the particular example of the exponential a-potential, v{a) = voexp{ga), 

(a) = 6 — ((/ + 6) [ 1 + Cl e ^, Ci = Cq {g + 6) , (63) 

which is positive for all ae3? if and only if Ci > 0, g < 0, g —6. Otherwise, it is not 
difficult to hnd solutions becoming negative on some intervals belonging to —cxo < a < +oo . 
In a moment we consider the case Ci < 0, g < 0 when the solution is rather complex. 

To simplify discussions of these complex solutions let us rewrite N/D formic 

=-gir - C 2 ){t + Ci)~^; C 2 = QCi/g, g = {g + 6), t = 6^°". (64) 

Then it is clear that for Ci > 0, g > 0 the nominator has one zero at r = —C2 and becomes 
negative for r > IC2I, where a > = ^“^ln|C2|. This means that is the branch point 

in the complex a-plane and there exists the second sheet of the Riemannian surface of the 
analytic function x(a) = having the ‘physical’ cut —00 < a < ao, with x(ao) = 0; 

on the upper edge of the cut x(a) is positive and on the lower edge it is negative. To hnd 
'ip{a) one should integrate x(a) along the cut; ao is the extremum of V’(a)- 

The singularities of x('^) corresponding to zeroes of the denominator are stronger because 
then x^(a) 00, but they are also integrable. The simplest such case is Ci < 0, > 0, 

when < 0 for r > |Ci| and is inhnite at r = |Ci| = where x^ ^ 9 (®i — «)~^- The 
same behavior has the special exact solution that can be obtained from Eq. fl63|l in the limit 
(7 = (6 + (7) —)■ 0: if we take |Ci| = 1 + ^ ao we hnd that the limit is (a) = 6 + (ao — a)“^ 
(this solution can also be derived from Eq. fl26|) for k = 0, I'{a) = —6). 

^®Note that the main properties of are defined by the signs of g and Ci, with the proper account of the 
restriction on the sign of (y^ — 6) according to equation (l28l) . 
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A configuration with more singularities emerges in the case Ci < 0, 5^ < 0. Then the 
‘positive support’ of the solution fl6T|) . where is positive, consists of two separate intervals: 

0 < r < |Ci|, ^2 < r < cxo, if ^ > 0; 0 < r < |C'2|, |Ci| < r < 00, if ^ < 0. 

One can see that there is a fundamental difference between solutions with different number 
of branch points, which also depends on their behavior (singular or regular). In fact, there 
exist hve types of the solutions: R (no b.p.), IR (1 regular b.p), IS (1 singular b.p.), 2SR 
(singular and regular b.p.). The support of the solutions IS and IR does not include large 
enough values of a and thus for them the metric has a hnite upper bound. We may call 
them gravitationally regular solutions. All other solutions do not have this property. 

Now it is possible to derive '0 (q^) by integrating along the ‘physical’ paths in 

the complex a-plane and thus to hnd the portraits of cosmologies with the simple exponential 
potential and k = 0. The same consideration can be applied to analyzing solutions for other 
potentials, e.g., v{a) = e^“P„(a). This class is rich enough to describe different scenarios, 
like inflation, bouncing, etc. 

4.2 Exact solutions for arbitrary v{a) and k 

Here we show that Eqs.([6])-([8]) can be transformed into linear equations for functions 'rf{a), 
^^(a), which can be integrated for arbitrary preassigned potential v{a) and with arbitrary 
curvature parameter k. It is more convenient to consider equations ([7j) and (jS]), which we 
transform to the a-picture using 0241) and other above dehnitions. Thus we derive two linear 
differential equations equations that can be solved with any given v{a): 

^ + 2(3 + c) 7f + v\a) = 0 , + rf + 2k = 0 . (65) 

da da 

Now, introducing new positive functions y{a) = rf and x{a) = we hnd for y{a) 

and x{a) the gauge-independent equations and constraint 

y'{a) + Qy{a) + v'{a) = 0, x'{a) + y{a) + 2fce“^“ = 0. (66) 

y{a) — 6 x{a) + v{a) + 6k = 0 . (67) 

Keeping in mind Eq. fl621) we write the solution of the hrst equation fl6B]) in the form, 

pa 

j/(a) = 6e“®"/(a) — v{a); I{a) = [Co+ / e^°‘v{a)], (68) 

J a— 

where Co and a- > —00 are real numbers. Then we hnd x{a) from constraint (1671) : 

x(a) = e-®"/(a) + (69) 

It can also be derived from fl6B]) but then the additional arbitrary constant must be hxed by 
constraint fE7)l . which shows that the integrals I (a) in fl65]) and fl69]) are identical. Now we 
can hnd the (gauge dependent) phase portrait of cosmology in the a-version, 

^ = a = e-“ [e-^“/(a) + k]^/\ (70) 

which gives t{a) by one integration over ‘physical’ cuts. As an exercise, one may derive 
explicit expressions for a{t) for the potentials v = e^°‘Pn{a). 
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The obtained solutions also give the exact expression for fh® case fc 7^ 0: 

X^{a) = y{c()/x{a) = [6e~^°'I{a) — v{a)][e~^°'I{a) + (71) 

which, of course, coincides with 062 p when k = 0. It is not difficult to check that this 
expression satisfies complete differential equation fl26l) if we take into account fl66l) . fl671) . 
Moreover, if we substitute into equation fl26P the expression for from ([69]), we 

obtain the well-defined equation for but it is much simpler to use the exact solutions. 

Anyway, the approach to analysis of the singularities and support of the most general solution 
is essentially the same as above - we first look for zeroes of the nominator and denominator 
and then find the cuts, on the edges of which we should integrate x('^) = 'ip'ioi) to find 
'll)(a). In principle, this result could allow us to establish a correspondence between the two 
versions. In practice, this is a not so simple task, which requires a careful investigation. 

At the moment, the best strategy is to compare physical results for interesting classes 
of potentials in the new version with those obtained for the well studied potentials of the 
standard -^-version. A more radical approach is to try to find, directly in the new ver¬ 
sion, potentials that describe physically interesting phenomena, like inflation, bouncing, or 
something else. As we demonstrated above, this approach can be significantly strength¬ 
ened by using in addition other inputs - various portraits motivated and supported both 
by theoretical intuitive ideas and observational data. The synergetic strategy of using the 
full mathematical structure outlined above - all gauges, versions, and inputs - looks like 
a promising global approach to isotropic cosmology that possibly could help us to return, 
sooner or later, in the higher-dimensional world of real physics. 

4.2.1 On replacing potential by kinetic energies 

It is usually supposed that the dynamical functions r] can be expressed in terms of the 
scalaron potential n('0), which is unknown and is usually chosen to satisfy some ‘reasonable’ 
properties providing a sort of inflation or other phenomena. It is supposed that the potential 
can be somehow derived in a future superstring theory or in present supergravity considera¬ 
tions. Our a-version may suggest a different approach to finding the potential - first guessing 
the scalaron kinetic ‘energy’ ri{a) and then deriving ^(0;) and .^(o!) as simple functionals of 
77. To derive the exact expression we simply integrate the relations fIBBD (forgetting that they 
are the differential equations for x{a), y{a)) and find: 

r°‘+ 

v{a) =—y{a) + QJ{a), x{a) = ke~‘^^ + J[a) ] J[a) = Ci + / 7/(a). (72) 

J OL 

Here we defined J{a) similar to I {a) introduced in (j68ll and with the same aim - to simplify 
considering the positivity conditions and asymptotic behavior for y{a), x{a), v{a). 

By the way, the expression for the fundamental function x(a) in terms of y{a), 

{<y) = y [ke~^‘^ + J {a)]~^ = —L'{a) + 1]~^, L{a) = \n J (a). (73) 

is simple and physically transparent. It demonstrates a deep connection between the scalaron 
kinetic energy and metric and permits to reconstruct all the portraits of our cosmology. The 
portrait a{a) is explicitly given by d72|) : the portrait 'ipi^a) = J x(a) can be derived from 
d73|) (for arbitrary k). As far as we can derive the inverse function alip) we also find V’(V’) 
and vl'ip) = h[a('0)]. Note also that y = —dL/dij) for fc = 0 to compare with Eqs. flT^ . flH^ . 
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One final remark on relation of this construction to inflationary models. The creators of 
inflation observed that very different potentials may give very similar inflationary scenarios. 
In fact, in the early Linde models the main idea was to hnd potentials that dehne inflationary 
‘motions’ of the inflaton. The basic ingredient was these specihc motions and the potential 
was an instrument to describe the model in a more standard held theoretic frame. Possibly, 
the expressions of the potential in terms of the kinetic energy of the scalaron depending on 
the metric may give a different, mathematically accurate realization of these ideas. 

Another option is to take as an input the main dynamical characteristic of the Universe, 
the Hubble parameter = x{a), together with x(a) determined by Eqs. fld^ . which is 

even simpler to use than Eq. fl73l) . This requires some work and time for customization but 
may give a new insight into the structure of cosmological models. One may call such an 
approach ‘constructive’ cosmology having in mind that the mathematical structure of the 
classical isotropic cosmology is an instrument transforming the input portraits or potentials 
into particular cosmological scenarios to be eventually confronted to the observational data. 
This idea looks not very appealing to theorists although some hints to ‘reconstructing’ the 
potential v{^jJ) from observational data, can be traced in literature. 0 

4.3 A fresh look at inflation and inflationary perturbation theory 

The standard approach to inflationary models does not contradict to the recent observational 
data, even if one uses various potentials and considers non-minimal inflaton coupling to 
gravity. However, it is not a completed theory as it actually describes the short period of the 
classical expansion and the emergence of quantum fluctuations on the classical background. 
At the moment, there exist a few other models also explaining main properties of the very 
early Universe but giving predictions that can in future be distinguished from those of 
inflationary modelsHj Our aim in this Section is to apply our analysis of global properties 
of the cosmological solutions and, especially, of the exact a-version solutions to identifying 
classical initial conditions and their consistency with the inflationary restrictions which we 
discuss in a moment. In this way, we hope to clarify the meaning of some specific features 
of the inflationary models and thus to make their comparison to other approaches easier. 

We begin with the standard conditions for inflation and the parameters accessible to 
measurements. In our language, the obvious necessary condition for inflation is |x('0)| 3> 1 
on a small interval of fj, or, equivalently, 1 on the corresponding large interval 

{ai < a < af), where «/ — ctj = W ~ 50 is the so-called number of e-folds (see, e.g., [3], 
[15]). Using equations flTT|) and fl66|l with fc = 0 we derive the exact relation^ 

- 2^/^^ = -x\a)/x{a) = y{a)/x{a) = x^{a) =2e = . (74) 

Here we dehne the new function e of one of the three variables (r, a, -0), which is supposed 
to be small, e <C 1, and approximately constant on the corresponding inhationary intervals 
{Ti,Tf), {ai,af), ('0/,'0i)- Under these conditions, we call it the hrst inhationary parameter 

After finishing this paper I became aware of Ref. m proposing to determine v{tjj) from some contempo¬ 
rary observational data on FLRW cosmology and, as a matter of fact, using for this purpose an a-version-type 
equation for the Hubble parameter. 

^®See, e.g., a discussion of some quantum aspects of the Universe creation in [48] that is related to a 
modernized version of the classical ekpyrotic-plus-bouncing dynamics. Our approach is applicable also to 
these more complex models. For a review of this an other alternatives to inflationary models see |49j . 

^®The standard inflationary models suppose that fc = 0, and in this Section we always keep this condition. 
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that is usually denoted by simple e. Applying now the exact constraint (1^ we derive 


f = i?lv(ii) = p{a)/v{a) = I (1 - i/3) ' - | • 


( 75 ) 


where the last approximation (important for the standard considerations) is valid if e <C 3 
6)0 The approximation can be applied only to positive potentials, and positivity is 
necessary at least in the inflationary domain. If > 1, the inflation stops and, for —)■ 6, 

the scalaron kinetic energy dominates over the potential one. If n > 0, all the solutions must 
satisfy the condition = l/y^ < 6; they can be only asymptotically close to the limiting 
curves y = ±\/6 (y = ±l/\/6 )• Near these attractors the parameter f becomes inhnite and 
the scalaron potential energy is negligible. As demonstrated in Section 3, this behavior of y 
is possible for large and small "0. In contrast, the condition |y| S> 1 is generally more difficult 
to satisfy in the ■^-version. It requires e: = — 1 in (IH)) and is possible for small enough ip, 
because \/6y(V’) ^ 2: = coth?/ can be large for small y only (see (H2]) . Here, we do not 

attempt to investigate this part of the '0-picture depending on unknown concrete potentials 
and initial conditions but describe an easier and more rigorous approach to inflation. 

Starting with the model-independent dehnition of inflation given above in frM)l - fl75|) we 
now show that the a-version gives the detailed picture of inflationary cosmology. In the 
standard approach to inflation, cosmologists use one more parameter fj dehned by 



y'{a) 

y{a) 


2(h - e), 


h-X^(a) 


y(a) “ di)~ ’ 


(76) 


where the second formula gives the exact expression of fj in terms of y(Q!) or y('0) = 1/y- 
The inflationary parameters can be derived in terms of v{a) by using solutions fl66l) - fl6^ and 
we can expand them in the derivatives h^"^^(a) applying the simple formula (/c = 0): 


N 


6x{a) = 6e-‘'“/(a) = Coc-^'^ + - (-6) 


—N^—6a 


g6a y(n+l) 


(a), (77) 


where Cq is an arbitrary integration constant dehning the unique solution [x{a] Co)^y{a] Co)], 
y{a) is given by fl66l) . The lower integration limit is chosen snpposing the effective potential 
v{a) vanishes in the extreme ‘quantnm’ limit a —)■ —cxo. A very useful property of the 
solution for the chaotic inflation is that the ratio of the hrst term to the rest is exponentially 
small. This means that in the ‘classical’ domain a > aj 3> 1 we may approximately use the 
solution with Cq = 0. Neglecting this term in the solution [x{a),y{a)], we hnd 


f (a) 


y{a) 

v{a) 




E(-i) 


6^v{a) 


(78) 


The properties of the series are dehned by the potential v{a) and, in principle, can be derived 
from the exact analytical solntion given by eqnations fl68l) . fl7T]) J^^I Thus we suppose the series 
dehnes a regular function y^ that is small in some interval 1 < ai < a < af. In addition, 
we assnme it is sufficient to take into acconnt just a few terms of the expansion, which give 

view of the factors (1/6)", which are most clearly seen in equation (1771) below, the inflationary 
conditions can actually be somewhat weaker, like e < 1/2. The exact bound depends on the potential. 

^^To avoid worries on the series convergence, we may have in mind polynomial potentials v{a) without 
zeroes in the domain of our interest, and thus the series has a finite number of terms. 
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the corrections to standard inflationary parameters. In this sense, Eq. fl78l) can be considered 
as a somewhat nnnsnal pertnrbative expansion. To relate it to the standard inflationary 
langnage we translate fl78l) into the -^-version by a perturbative construction of y. 

The lowest approximation for y is given by the first term of the snm in ([78]): 


v\a) 

v{a) 


-I'(a) = -X I'W = -X , 


x = xi = -I’W ; hi = 


v"W 


(79) 


where iji is derived from fl76|) with x approximated by Xi and, in the same approximation, 
2e = Xi = (^0^- Here we neglect the terms with n > 2, which can be calcnlated recnrsively 
and are sums of the monomials nr=i(^^*^) with ^ i/cj = n, as can be seen from relations: 


-,(«+!) 


v[a] 




da via) 


v[a] 


'^ = l'(a), 'X = [«+ (iX,...; (80) 

Via) V 


we need not write the same expressions for v^'^'^i'ip)/vi'ip) in terms of . It must 

be emphasized that the ‘logarithmic’ representations for x^ and x (invariant under v-scaling, 
i.e., V I—)■ cqv) are strictly valid only for solutions with Co = 0 in fl77|l l^ We have v-scale 
invariance if and only if Co = k = 0. This is a most general definition of inflationary x ■ 
Now let us And a few first terms for the expansion of x, which give corrections to the 
standard inflationary parameters. Defining da = d/da, dp = d/df) and recalling that 
da = xdp , we express the operators d” in terms of powers of dp and of x: 


da = x[xd% + (dpx)dp] , dl = x[x^d\ + 3x(dpx) dj, + dp(xdpx) dp],... d” = xDn , (81) 


where the operator x^n can be calculated recursively. This operator is homogenous and of 
the same order n in both the variables dp and x- This means that the n-th term of the 
perturbative expansion (1781) is of the order 2n in dp, which is also of order one, (179jl . 0 By 
substituting these operators into x^ from (j78|l we And, after division by x, the exact relation 




Dn * vj'i/) 
^ (—6)"“^n V 


D2*V 
V l'{4>) 


j[l'' + [I'Y] + (dpx), 


(82) 


where * n('0) denotes the action of the differential operator on vif)), which is illus¬ 
trated by the D 2 example (note that x/i' = — 1 + 0(Xi) ^'Cid thus I' in denominators is not 
dangerous). It is now easy to derive the second approximation X2 by replacing in fl82|) x by 
Xi and neglecting the terms of the order n> 3: 


X 2 = -I'if’) 




1 + -121" + HX]} = -I'X) 


1 + f + 0{xt) 


(83) 


To calculate the third-order term we substitute X2 for x ici ll8^ and take in the sum the 
terms up to n = 3. This may be a simple exercise for a reader, o 

To illustrate perturbation theory we write the corrections to the inflationary parameters, 
which follow from (18^ for x and x = 1/x- The corrected ei is 62 = C[1 + 2/"/3]; for the 

^^The general solution has the non-scaling non-logarithmic term Cq e“®“/u(a) in the (1751) . 

^^We now see that the cosmological ‘parameters’ e and fj are of the second order in dp while the funda¬ 
mental first-order parameter of our perturbation theory must naturally be xi = I'. 

simpler exercise is to show that the condition of linearity = 0, with x replaced by xi, defines 

bilinear potentials u(■)/’) = uq ('b + 'bo)^- 
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simplest potentials v = Vo'ip"^^ we find in this approximation e = ei[l — 2ei/3N + 0(e^)], 
where ei = Xi/2 = 2N‘^ and |2ei/3A^| must be small. Not much more difficult is to hnd 
172- From fl76|l . fl791l80|) we see that 171 = V + (Z')^ and, using ([83]), we hnd 

m = xl- X2 = fii + (1/3) [ 2 Z(2) + /(i)/(3) _ ('7(2))2| _ 

Taking for illustration the potential v = vo'ip'^^ we easily hnd that 

r)i = 2N{2N - 1) V’-' , r)2 = hi{ 1 - hi [ (4iV - 3)/ 3 {2N - 1)' ] } , 


which is h2 = 'i7i( 1 — hi/3) for = 1. With X — X 2 , fhe number of e-folds is 


Ctpf 


Ne = 


d^PxW = 


fpi 


hpf 


d'lp 

W)[ 


l-f + 0{x\) 




f 

i ’ 


( 86 ) 


where the hrst term is the standard result, and W = [i/^/4iV + Ini^/d]/ when v = vo'ip'^^. 
Now, using Eq. flS^ . we can approximate equation in the inhationary domain S> 1, 

dz/dx = {z^ — l)[zu{x) + 1] = — 1) [u'{x) + ...]., (86) 


see notation in (H2|) . Then, recalling fHi|) - (H5|) and taking z{x) = — coth(7/), we hnd that 


dy/dx = 1 — u{x) coth{y) = u'{x) + ...; z = coth(xi/2) = — cothM(a:), (87) 


where we have chosen an arbitrary integration constant zero (to have z ^ 00 when m 0). 
Of course, this relation does not tell us anything new: returning to y and ijj we simply hnd our 
hrst-order relation y = —\/6 tanh(/'/\/6) = —V(7p). However, this observation may signal 
that a perturbation series for the inhationary solutions can somehow be derived directly in 
the "0-version, most probably, for special solutions like fl55|) . (156|) . It is instructive to look at 
a simplest asymptotic expansion of y{x) in powers of 1/a: ~ ^(a;) ~ for n = vqx"^^: 

00 

y{x) = ■, yi = N, 7/3 =, (88) 

0 

which is obviously similar to inhationary perturbation expansions. This approach can prob¬ 
ably be generalized to potentials for which u{x) is expandable in powers of 1/x. Similar 
asymptotic expansions can be used as approximations of y{x) near the hnite-range singu¬ 
larities like u{x) ~ l/(x — Xi) produced by the poles of the Higgs-type potentials v{x). For 
inhationary potentials depending on powers of for which u{x) = Un , we can 
analogously exploit the expansions y = . 

Note that the hrst-order terms of our approximations to the exact inflationary a-solutions 
reproduces the essentially approximate formulae used in the inflationary models, while our 
higher-order terms give signihcant corrections to them. We also stress that all three inha¬ 
tionary parameters can be expressed in terms of one fundamental function = 2e which 
describes all properties of the model. The second parameter, f), can be explicitly expressed 
in terms of by Eq. fl76|) . which shows that fj must be small as far as y and d^x small. 
The parameter W can be expressed as an integral of y = y“^ over an inhationary -(/-interval 
and also is produced by y. Of course, in the present-day observations one measures not 
the y-functions but some averages of e, i), which nevertheless provide us with an important 
information on physics of inhationary models encoded in y function. 

Note hnally that the a-version of inhation gives a direct and smooth transition from the 
‘quantum’ domain a < 0 to the classical one, where a 3> 1. This picture looks surprisingly 
simple and coherent although further work is needed to clarify its relation to the standard 
cosmology. Even more fascinating might be attempts to leak out of the classical to quantum 
domain using ideas as simple as the classical dynamical cosmology. 
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5 A short summary 

Here we briefly summarize the main results and problems that should be clarihed in future. 

1. The cosmological dynamical equations are formulated in different gauges and versions. 
We illustrate relations between them on simple solutions and by integrable models. 

2 . The general properties of gauge independent x-^Quations fl26l) - flTr|) . describing the main 
{a,if) portraits of isotropic cosmologies, are established in a and if versions. 

3 . Equations fl2^ - fl30|l allow us to derive the complete solution if x(a) or x{if) are known. 
Taking into account equations (l3Tll we can in addition derive v{a) or v{if). 

4. We discussed different ways to determine cosmologies not using potentials. A most natural 

one seems to hrst derive using (l32|) . with the Hubble function as an input. 

5 . Although the x-equations depend only on v'{if)/v{if) and are thus insensitive to the sign 
of v{if) = v{a), this sign is critically important for global properties of the solutions. From 
fl28p . (HID . (ITS]) it follows that the solutions in the intervals with v{if) > 0 are isolated from 
those in the intervals with v{if) < 0 and must be studied separately. 

6 . We mostly considered potentials not changing the sign and studied in detail models 
with positive potentials for which inflationary scenarios are natural. We also can use and 
actually used our solutions and their expansions near the points ifo where v{ifo) = 0 and 
thus vfif)/v{if) behaves as {if — ifo)~^ —t ±cxo. This is a problem in the -^-version because 
{xf' — 6) may change the sign with the potential, as follows from fHT|) . fl75D . But in the 
a-version it is no problem at all, as can be seen from from expression flB2]) for 

7. Probably, the most important results are presented in Section 4, where we have found the 
exact solution of all equations for arbitrary v{a) and k. The necessary condition for inflation 
is x^{<^) < 6 {Qx^iif > !)• To derive from x(a) standard inflationary scenarios we hrst 
suppose that the spatial curvature vanishes, k = 0. Then, by hxing the arbitrary integration 
constant, Cq = 0, we preserve the n-scale invariance of inhationary solution y and derive its 
expansion from Eq. flS^ as a sum, the n-th term of which for n > 1 has the form: 

2n 

-l'{if)'^Cn{ki, ...,k2n)Y\_U^'^W]’^'y where ^ ifci = 2n , ki>0. 

k i=l i 

This inflationary perturbation expansion can be obtained by the well-dehned recursive alge¬ 
braic iterations and gives higher-order corrections to the inhationary parameters e, rj ,Ng . 

8. When v{if) < 0 and thus x^(q;) > 6, 6x^{if) < 1, it is also convenient to use expansions 
of x{if) when it is small or close to 1/6. In the last case we have derived asymptotic ap¬ 
proximation fl5T|) for -0 —>• cxo valid for a broad class of potentials v{if). We should mention 
an interesting one-parameter class of ‘bouncing’ solutions (see fl59D and fl60D ). which exist 
when v'{if)/v{if) ~ 1/V’; cind a special solution fISHD that probably is a separatrix. The 
global picture of solutions with such properties are of great interest for ekpyrotic-bouncing 
scenarios and must be studied in future. 

I believe that visualization of these structures, drawing the {a, if) portraits, and using per¬ 
turbative expansions for concrete inhationary, ekpyrotic, bouncing and other, more strange 
isotropic cosmologies may stimulate their better theoretical understanding. 
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6 Appendices 

6.1 On general isotropic cosmologies 

A fairly general dimensional reduction of the Einstein gravity coupled to a scalar field ifj, 
which gives all possible spherically symmetric cosmologies is described in [23], [36]. Following 
this procedure we derive the effective two-dimensional Lagrangian describing spherical static 
states, cosmologies, and waves. The starting point is the two-dimensional metric of the 
spherically symmetric space-time (we usually denote = ip and call it the dilaton field), 

dsl = e^°‘dr‘^ + 0) — e^'^dt^ -|- 2e^^drdt , (89) 

where a, jd, 7 , 6 depend on t, r and dVl?{0^ cf)) is the metric on the 2 -dimensional sphere 
Then the two-dimensional reduction of the four-dimensional Einstein gravity coupled to a 
scalar field 'ip is well known (here the prime denotes differentiations in r and the dot - in t): 

^( 2 ) _ ga+2/3-7^^2_2^2_4^^^_g-a+2/3+7^^/2_2^/2_4^Y^_ga+2/3+7p.(^)^2fce“+^, (90) 

where fc = 0, ±1 is the standard curvature parameter, which vanishes for spatially flat 
cosmologies. This Lagrangian is obtained from the complete four-dimensional one, ([I]), by 
omitting the last term in metric fl89p and extracting the total derivative terms in 0901) : 

A/:(2) = _2[(e^)'e2/^-“ + 

(the omitted time derivative part is derived by replacement dr -H- idt and a y-)- 7 ). 

Let us also recall that fl90|l is a gauge theory with two constraints: the total energy 
and momentum vanish according to the equations of motion derived by variations in all the 
variables. The origin of the constraints can be related to independence of the Lagrangian of 
the derivatives d and 7 '. It follows that an^l become Lagrangian multipliers in the 
one-dimensional static and cosmological reductions, respectively. A more rigorous treatment 
requires applying the ADM Hamiltonian formulation, [50]. 

Variations of this Lagrangian give all the equations of motionl^ except one constraint, 

- - ^fd'+ dl3'+ (91) 

which should be derived before we omit the 5-term in the metric (taking the limit 6 —)■ —cxd). 
All other equations of motion can be obtained from the effective Lagrangian fl90|) . 

Now, the distinction between static and cosmological solutions is in the dependence of 
their ‘matter’ field pj on the space-time coordinates. We call static the solution for which 
pj = 'ip{r). If p) = 'ip{t) we call the solution cosmological. There also exist the wave-like 
solutions for which a, /3 ,7 and pj may depend on linear combinations of t and r but we here 
do not discuss this possibility. For both static and cosmological solutions the gravitational 
variables in general depend on t and r. This is important for the embedding the solution 
into higher dimensional theory but here we may forget about the dependence of cosmological 
solutions on the space coordinate. 

^®The equations of motion are the standard Einstein equation in the spherical coordinates. We need not 
write any Lagrangians to apply to them further separations of variables. However, introducing effective 
Lagrangians and Hamiltonians is extremely convenient, even in the classical environment, and will become 
indispensable if we turn to quantizing them (e.g., m)- 
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To obtain the one-dimensional eqnations we make farther redactions by separating t and 
r. It is clear that to separate the variables r and t in the metric we should require that 

a = ao{t) + ai{r), /3 =/3o{t) +7 = 7o(t) + 7i(^), (92) 

Inserting this into the equations of motion one can hnd the restrictions on the gravitational 
(and, possibly on the matter) variables that must be fulhlled. The details can be found in 
[23] , where one can hnd the complete list of the static and cosmological spherically symmetric 
solutions. The naive cosmological reduction (that supposes all the helds to be independent 
of r) does not give the standard FLRW cosmology with one scalar. As was shown in [23] 
(see also the earlier paper 1201). one of the possible systems of conditions for separating the 
variables in the Einstein equations or in the Lagrangian is the following 

a = $, 7 ' = 0 , /?i" k = 0 , 2/?i" + 3/3[^ - k = 3k , (93) 

where the hrst two follow from Eq.(j92|). The constant k in the equations for fii is propor¬ 
tional to the 3-curvature of the space-time, and the third equation in (l93ll is the isotropy 
condition. Any homogeneous isotropic cosmology must satisfy all four conditions^ Both 
the homogeneity and isotropy conditions follow from one equation 

(94)/3(^ - = fc, (94) 

which can easily be solved for all values of the parameters. Using Eqs. (l93|l . (IMll we get the 
standard effective Lagrangian (|2|). We see that for naive reductions the isotropy conditions 
in (|93|) can be satished only if = 0 and that the hrst condition is not dictated by (l9l|) . 
Therefore, naive reductions give, in general, homogeneous non-isotropic cosmologies. 

In case of fc = fc = 0 we can show that in the arbitrary naive cosmology with /c = 0 there 
exists a class of ‘isotropic’ solutions satisfying the condition a{t) = a{t) — /3{t) = const. It 
is easiest to demonstrate this in the Hamiltonian gauge 7 = a -|- 3/3 in which: 

a = a — /3 = k , 2ci = n(V’) 2'^ = —. (95) 

When /c = 0, there exists the integral of motion & = Cq and thus (3 = a — C^t. The solutions 
belonging to the class with Uq = 0 are isotropic. Moreover, if the potential is exponential, 
i.e. v'kip) = gv^ip) there appears one more integral, xj) + ga = Ui, and the model becomes 
integrable. Then the equation for a{t) can be easily reduced to the Liouville equation and 
explicitly solved. The result can be presented in a gauge independent form xjj = /(q;)@ 
This model of a relation between isotropic and anisotropic cosmologies is, of course, 
unrealistic. Whether some anisotropic cosmologies may have physically interesting isotropic 
limits is an interesting problem which is not discussed in this paper. With an additional 
scalar held, an evolution to isotropy looks possible, but then one should consider the complete 
system of three equations discussed in [36]. This problem requires a separate investigation. 

6.2 Integrable example of equation for xO^) 

Here we hnd the potential for which the k = 0 reduction of Eq. fl64ll can be exactly solved 
and derive the correspondent solution, the general structure of which turn out similar to that 

^®Note that we here neglect inessential constant factors and have chosen ai = 71 = 0. 

^^Using the two integrals we find the equation for a looking like d = exp[o(t) + bt]. Denoting a + bt = x{t), 
we find the Liouville equation x = exp a:. This allows to find t{a) and exclude t from the second integral. 



of the general a-solution. We first find a special solution Za{x) of (H2]) by supposing that 
u{x) Za(x) = a — 1 where a will be determined later and Za satisfy z'^ = a{zl — 1). Obviously, 
Za = — coth^(aa;) where e = ±1 and thus u{x) = (1 — a) tanh^ ax, where x = x — Xq. Putting 


from now on xq = 0 we hnd the corresponding potential v{x) (recall ( 02 ])), 

h(x) = vo{ I |/2 ) , £ = d=l. (96) 

By the way, u{x) = v'{x)/2v{x) satishes the simple equation 

u'{x) = a(a — — (a — 1)^]. (97) 

Now, let us try to hnd the general solution by the substitution 

z = y + Za= y + {a — l)/u = y — cotlP(ax). (98) 

It follows that y{x) satisfy the Abel equation, 

y' = uy^ + {a + 2u Za) y"^ + {2 a Za + u zl - u) y, (99) 

which can be explicitly integrated if a + 2uza = 0] this is possible if a = 2/3. Then, 

y'= uy^ — {u + l/3u)y, 1/3m = coth^(2x/3) (100) 

and w = y~‘^ satishes the linear equation: 

w\x) = 2{u + l/3u)w-2u. (101) 

We write the general solution for e = —1 which is regular at x = 0 (a = 2/3), 

w{x) = cosh^(ax) [cosh( 2 ax) + Co sinh( 2 ax) ], z{x) = {w{x))~^ + tanh(ax). ( 102 ) 


The solution for e = +1 will be regular if we choose the minus sign of the root in fll02p 
to cancel the singularity of Za = 1 / 3 m ( x ), see fl96]l . The function w{x) is positive for all x 
ifC^ < 1 and is equal to cosh^(2x/3) exp(±4x/3) for Cq = ±1. Similarly to the general 
solutions in the a-picture, we have here the positivity problem and the second sheet also 
emerges if Cq < —1. In any case, by integrating z{x) we can derive the explicit expression for 
q;('^) that, in principle can be inverted to obtain '0(a). Thus '0(a) and 'y(a) = 'y['0(a)] will 
implicitly depend on two arbitrary parameters - one is Cq and the second is the integration 
constant ' 0 o emerging in deriving 0 (a). 

6.3 On versus v{a) in simple solutions 

Here we consider important examples of deriving the potential v both in a and 0 pictures. 
We discuss in some detail the problem of hnding potentials for which there exist some simple 
solutions of equations (1T0|1 - (1T^ . where we omit bars and denote r-derivatives by dots: 

2i + ri‘^ + 2ke-^^ = 0, 2fi + Gr]^ + v'{^p) = 0 , (103) 

V{%[)) = — rf — Gke~‘^°'. (104) 

The simplest and physically interesting solutions constructed on one of two ad hoc guesses: 
1. ^ = Co OT 2. T] = Cq. The key idea is to solve the potential-independent equation, to 
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derive x('^) or x(V’) giving a('^), and then to use constraint 01041) for finding n(V’)- The 
aim is twofold: hrst to learn something on a rather nontrivial relation between the standard 
and a ‘versions’ and, second, to hnd simplest ‘dual’ potentials ^(-0) and v{o). The result of 
these simplest, almost trivial considerations looks somewhat unexpected, especially, in case 
of non-vanishing curvature term 

We consider in some detail only the hrst, simpler solutions. When Co = 0 we easily hnd 
d(r) = ^ = ,^0 , a(r) = - tq) , r] = ko = -2k). (105) 

The partial map x(a;) and the corresponding ijj{a) are very simple in this casein 

X{a) = ^ = I = ^ ^ - ^o) = y x{a) = e"“ ; xW = -i . (106) 

It follows that the potential in a and in 'ijj versions is very simple, 

v{a) = 6 + 2 kl = v [a{i))] = 6 + 2 “ V’o)^ = v{'ip). (107) 

It is not difficult to check that all dynamical equation are satished in all versions. Note that 
the arbitrary parameters have diherent meaning - ,^o is the Hubble constant while integration 
constant 'ipo is unimportant as V’ is in fact dehned up a shift. 

The case ^ = Co = —Vo/2 is solved quite similarly and we skip details. We have: 

= r]l + kl , ^ = Co(r-ro), a = Cq {t - tq)'^/ 2 + ao , ^ = -Vo V^o - « • (108) 

Now it is easy to derive the portrait and the simple expression for its limit for ko —)■ 0, 

= y x{oi) = J da {1 + klvo"^ {oio — 2y/ao — a , (109) 

and easily hnd v{a) and for the last we only write fco 0 limit: 

v{a) = Qvo («o -a)-vl + ‘2kl = v[a{^/J)] (3/2) Vo O’ - V’o)^ - Vo ■ (HO) 

This limit gives the same potential ~ obtained for the ^ = 0 example and this also is true 
for the V = 0 case. The condition k = 0 signihcantly simplihes solutions, especially, in the 
a-version (see Section 4). Equations fllOOp - flllOp with ko 0 and similar ones for the case 
V = Co ^ 0 give more complex potentials. Though the problems are of physics interest we 
cannot them in the present paper. 

It would be of interest to hnd and study potentials giving physically interesting solutions 
and at the same time analytically accessible in the a-picture. Such an approach looks 
viable and deserving careful elaboration but it will require much deeper understanding of 
the analytical structure and the precise meaning of 'i/>{a) as well as of the a-version as a 
whole. For these reasons, we are compelled to leave this task to future investigations. At the 
moment this problem is solved in the frame of inhationary perturbation theory. Section 4. 
As a next step in this direction, one may analyze from this point of view the examples of this 
Appendix and explicitly integrable models briehy discussed in the main text of this paper. 


^®In fact, we find only the part of the complete portrait that corresponds to the chosen solution. As is 
demonstrated in the main text, we can reconstruct the complete portrait for any potential v{a) in the a 
picture but even in this simplest example we do not know the portrait for the corresponding potential vi^jj). 
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